Bound states in Andreev billiards with soft walls 
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The energy spectrum and the eigenstates of a rectangular quantum dot containing soft potential 
walls in contact with a superconductor are calculated by solving the Bogoliubov-de Gennes (BdG) 
equation. We compare the quantum mechanical solutions with a semiclassical analysis using a 
Bohr-Sommerfeld (BS) quantization of periodic orbits. We propose a simple extension of the BS 
approximation which is well suited to describe Andreev billiards with parabolic potential walls. 
The underlying classical periodic electron-hole orbits are directly identified in terms of "scar" like 
features engraved in the quantum wavefunctions of Andreev states determined here for the first 
time. 
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I. INTRODUCTION 

The interface between a normal-conducting (N), bal- 
listic quantum dot and a superconductor (S) gives rise 
to the coherent scattering of electrons into holes, 1 This 
phenomenon, which is of great experimental 2 *^ and the- 
oretica6SiL2i2ii2iiiii2ii^ interest, is generally known as 
Andreev reflection™ A N-S hybrid structure consisting 
of a superconductor attached to a normal cavity is com- 
monly called an Andreev billiard. For a recent and com- 
prehensive review of these systems see, e.g., Ref. HEI and 
references therein. 

An important quantity for Andreev billiards is the den- 
sity of states (DOS) which has been studied theoreti- 
cally by many authorsi 16 i 17 i 18 i 19 i 20 i 21 i 22 i 23 i 24 Its behavior 
close to the Fermi energy depends on the classical dy- 
namics found in the isolated normal conducting cavity: 
In the integrable case, the DOS is proportional to the 
energy, while for a classically chaotic normal cavity, a 
minigap (which is smaller than the bulk superconductor 
gap Ao) develops. Moreover, depending on the geometry 
of the normal cavity one can observe singularities in the 
DOS. A well-known example of such a behavior of the 
DOS is an Andreev billiard formed from a normal metal 
film attached to a superconductor studied long ago by 
de Gennes and Saint- James™ Recently, similar singular 
features of the DOS have been found in other Andreev 
billiards. 20 i 21 i 22 i 24 i 26 

To calculate the DOS for Andreev billiards many re- 
searchers have successfully applied the Bohr-Sommerfeld 
(BS) approximation, 15 i 16 i 17 i 18 i 19 i 20 i 21 i 22 i 24 i 26 i 27 i 28 It was 
shown that the DOS can be related to the purely 
geometry-dependent path length distribution P(l) which 
is the classical probability that an electron entering the N 
region at the N-S interface returns to the interface after 
a path length I. 



In the present work, we discuss on the one hand quan- 
tum and semiclassical results for Andreev billiards with 
walls which are mediating a smooth transition between 
the interior of the cavity and the region outside. On the 
other hand, we calculate, for the first time, the wave- 
functions of eigenstates in such Andreev billiards. As a 
prototypical example for an Andreev billiard confined by 
soft walls, we choose a parabolic wall profile, as shown 
in Fig. ^ Such soft-walled Andreev billiards are of in- 
terest because they describe the potential profile found 
in studies of typical quantum dots realized by remote 
surface gates. 29 Soft walls are a quite realistic approxi- 
mation which extends previous theoretical work where 
infinitely high walls were employed. We demonstrate 
that with suitable adaptions, the BS approximation can 
describe N-S hybrid systems in which the electrostatic 
potential in the normal conducting cavity is modelled 
by a parabolic potential. Recently, Silvestrov et al™ 
have presented a quasiclassical study of Andreev billiards 
containing smoothly varying potentials inside the normal 
dot. However, no exact quantum mechanical calculation 
of the energy levels of such Andreev billiards has been 
performed. We succesfully tested the predictions for the 
DOS obtained from our BS approximation by compar- 
ing them to results found from exact quantum mechan- 
ical calculations using the Bogoliubov-de Gennes (BdG) 
equation. 

A better understanding of the BS approximation can 
be gained by the analysis of electron and hole wavefunc- 
tions of Andreev states which we determine by the mod- 
ular recursive Green's function method (MRGM)™^ 
Based on recent theoretical studies^ 3,34 the wave function 
patterns observed in Andreev billiards may be scanned 
by measuring the tunneling conductance of such systems. 
We show that the wave functions feature enhanced den- 
sity over continuous families of classical periodic orbits 
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FIG. 1: Andreev billiard with parabolic walls at (a) the wall 
opposite to the lead and (b) at all sides except the one the 
superconducting lead is attached to. The superconducting 
area is shaded. 



( "bundles"—) . These bundles of classical orbits give rise 
to peaks in the pathlength distribution P(l) and, as a 
consequence, the DOS shows singular behavior at cer- 
tain energies. Moreover, a close similarity between the 
electron and the hole components of the eigenstates is 
observed in line with the semiclassical picture of the An- 
dreev retroreflection process. 

This paper is organized as follows. In Sec. [D] the 
method for calculating the quantum mechanical eigenen- 
ergies and eigenstates is presented along with our ex- 
tended semiclassical approach. In Sec. IHII our numerical 
results are presented. Finally, the paper concludes with 
a short summary in Sec. IIVI 



II. METHODS 

In this section we outline our quantum mechanical 
method for determining the energy levels and the cor- 
responding eigenstates of the N-S hybrid systems shown 
in Figs. [Ha) and (b). Moreover, we present the semiclas- 
sical approach for the density of states. 

The normal region of the N-S hybrid system is a bal- 
listic, rectangular cavity. The length of the side of the 
cavity parallel to the N-S interface is L. For the spatial 
dependence of the confining potential V(x), we consider 
in this work two cases: i) one parabolic- type soft wall is 
placed opposite to the N-S interface and all other sides 
of the rectangle (apart from the N-S interface) are hard 
walls (see Fig. QJa)), i.e. 



V(x) 



aE F (x — xq) 2 & (x — xo) if 0<y<L 
oo elsewhere 



(la) 

and ii) three sides of the normal cavity are confined by 
a parabolic- type soft wall (see Fig.QJ&)) and hence V(x) 
is given by 

V(x) = aEp (x — xq) 2 Q(x — xo) + y 2 S(— y) 

+ (y-L) 2 0(y-L)], (lb) 

where in both cases a is a parameter controlling the 
steepness of the soft wall, O is the Heaviside func- 
tion, while the parameters xq and L fix the position of 



the potential. The superconducting region of width W 
(W < L) is attached to the remaining side of the cavity. 
An ideal interface between the normal and the supercon- 
ducting region is assumed, i.e. the effective masses and 
the Fermi energies are the same in the N and S regions, 
and no tunnel barrier is present at the interface. 

The superconducting pairing potential Ao is constant 
in the S region and zero in the N region. This approxi- 
mation is valid if the superconducting coherence length 
£ is small compared to the characteristic size of the rect- 
angular cavity. 8 Thus, the pairing potential is given by a 
step function model A(x) = |Ao| 0(#ns — x )i where the 
coordinate system x = (x,y) is chosen such that the N-S 
interface is located at x^s — 0. 



A. Quantum mechanical solution 

The quantum mechanical description of the N-S hybrid 
system is given by the BdG equation^ 



A* 



-m 



v(x) 



v(x) 



(2) 



where Ho = p 2 /(2m*) + V(x) — is the single-particle 
Hamiltonian with the confining potential V(x) that de- 
fines the normal dot and m* is the effective mass. The 
electron and the hole part of the quasiparticle wave func- 
tion are denoted by u{x) and w(x), respectively, whereas 
e is the excitation energy of the quasiparticle measured 
from the Fermi energy Ep. In this work we study the 
bound states of Andreev billiards, i.e. the eigenenergies 
s are in the range < e < A. 

The energy levels S{ of the hybrid system are found 
by matching the wave functions at the N-S interface. To 
construct the wave functions satisfying the Bogoliubov- 
de Gennes equation in the two regions we follow the 
methods of Beenakkeri^i and Cserti et al.. 21 The wave 
function in the normal region can be expressed in terms 
of the scattering matrix S(e) of the open system in 
which the superconductor is replaced by a normal lead. 
This scattering matrix is calculated by the modular re- 
cursive Green's function method developed by Rotter 
et al. 31,32 Note that this method to study Andreev bil- 
liards is, within the model assumption outlined above, 
exact. In particular, it does not rely on the usual An- 
dreev approximation^ i.e., Ao <C £^f, and the assump- 
tion of quasi-particles whose angle of incidence/reflection 
are approximately perpendicular to the N-S interfaced 

The integrated density of states (in the following called 
state counting function) can be obtained from the energy 
levels Sj as 



(3) 



In our numerical calculations discussed below Nqm(s) 
can be obtained directly. Therefore, it is straightforward 
to compare the state counting functions obtained from 
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the exact quantum mechanical calculations with that 
from the BS approximation. The numerical differenti- 
ation to determine the DOS can thus be by-passed. 

We now sketch the method for calculating the wave 
functions in the two regions. From the matching condi- 
tions one can find the expansion coefficients b e n , b% of 
the wave functions in the S region, and c^ ,=t , in the 
N region in terms of the right (+) and left (— ) moving 
plane waves given by 



Semiclassical treatment 



Xn(x;£) 



2m* 



+ixk x 



hk x . n W 



n sm(yk y , n ) : (4) 



where k y , n = nir/W and k 2 n = 2m*E/h 2 



ky n are the 



transverse and longitudinal wavenumbers respectively, 
and W is the width of the N-S interface (here the wave 
functions are flux normalized) . The wave function in the 
S region is given by 

Us (x) = [lKXn(^E + )+Yb h n xt(^E-)] (5a) 

n 

v s (x) = ^[^Xn(x;^ + )+^x+(x;^-)], (5b) 



where 7 = A/(e + W /S 2 -s 2 ), E ± = E F ± WA 2 -e 2 , 
and the asterisk denotes complex conjugation. Note that 
to ensure the boundary conditions at x —> —00, only the 
right (left) moving plane waves x+(x;£?~) (x^(x;£? + )) 
enter in the expansion of the wave function in the S re- 
gion. 

Calculation of the wave function in the normal dot 
requires two steps. At first we obtain the wave function 
at the N-S interface using the expansion coefficients 
and cjj'^ calculated from the matching conditions at the 
N-S interface. Secondly, the wave function inside the dot 
can be written as 



U N (*,S) = ^<' + ^n(x,£), 
n 

v N (x,e) = ^c£~Vn(x,-e), 



(6a) 
(6b) 



where the ^ n (x, e) are the scattering wave functions of a 
particle inside the dot which can be obtained by project- 
ing the retarded Green function G(x, x', E) of the cavity 
onto the incoming wave^ 



^n(x, e) 



dy' G(x,x',E F 

Jo 



where rj n (y) y/2/W sin (^) and v x , n hk x , n /m*. 
While the wavefunction in the superconductor is given 
as a sum of analytically determined functions^ in the 
continuum limit, the wavefunction in the normal con- 
ductor is determined numerically on a tight-binding grid. 
In spite of these two very different approaches, we were 
able to fulfill the matching conditions with remarkable 
accuracy. The latter is a measure for the degree of con- 
vergence of the tight-binding grid calculation towards the 
continuum limit. 




FIG. 2: A rectangular Andreev billiard with one parabolic 
wall defined by Eq. (|lafl (darker shading in the N region cor- 
responds to higher potential). One typical Andreev orbit, 
consisting of an electron part (dashed line) and a hole part 
(dotted line), is shown. 

Over the last decade, the Bohr-Sommerfeld 
approximation for the smoothed density of 
states has been successfully applied to Andreev 
billiards, 15 ! 16 ! 17 ! 18 ! 19 ! 20 ! 21 ! 22 ! 24 ! 26 ! 27 ! 28 In the case of 
a normal dot confined by one N-S interface and in- 
finitely high potential walls, the integrated density of 
states or smoothed state counting function in the BS 
approximation reads 



^BS(^) 



00 



P(l)dl, 



(8) 



where the integer part of M = k-pW/ir is the number of 
propagating modes in a lead of width W at the Fermi 
energy E? = h 2 k^/ '(2m*) (here fcp is the Fermi wave 
number) and P(l) denotes the path length distribution, 
i.e. the classical probability density for electrons entering 
the billiard at the N-S contact to exit after a path length 
I. Finally, l n (e) is given by 



U7T + arccos 



£ 



(9) 



with vy the Fermi velocity. In order to derive Eq. Q, one 
assumes that the hole retraces the path of the electron 
(retracing approximation, see Fig.EJ). This allows one to 
quantize the periodic orbits created by two subsequent 
Andreev reflections using the BS quantization condition 
for N-S systems: 



(Pe - Ph)dq : 



27rft (n + t 
4 



>)• (10) 



In Eq. (ITOl) a stands for the Maslov index (for details see, 
e.g., Ref. 39), 4> = - arccos(6:/A) represents the energy 
dependent phase shift resulting from Andreev reflection 
at the N-S interface and T is an arbitrary path of geo- 
metric length l g connecting one point at the N-S interface 
with another (see e.g. the one indicated in Fig. EJ). The 
integral in Eq. ([!()[) for a normal dot with hard walls re- 
sults in S e — Sh = 2eljvY from which Eq. (0) follows. 

We now extend the BS quantization to soft walls con- 
structed from harmonic potentials. Due to the form of 



the confining potential V(x) = V(x) + denned in 
Eq. (JU, the component of the momentum parallel to the 
soft wall is a constant of motion. Therefore the action 
(Eq. (O) for the part of T which lies in the V(x) ^ 
region can be decomposed as S e — Sh = S\\ + where 
5|| (Sj_) involves the momentum components parallel 
(Pe,|| — Ph,\\) an d perpendicular (p e ,± ~Ph,±) to the wall. 
Calculation of S\\ is trivial since p ej y — p^j| is constant, 
resulting in 5y = (p e — Ph)l\\ sin ^ where $ denotes the an- 
gle between p e and p e n or between p^ and p^n, while h 
corresponds to the displacement parallel to the soft wall 
which is approached by the particle. S± can be evaluated 
for the parabolic potential (Eq. JJJ) as 



S±=1T f!A- e cos 2 d. (11) 
V aEp 

Note that within the retracing approximation the vec- 
tors p e and ph at a given position x are antiparallel to 
each other while their magnitude is different for e ^ 0, 
i.e. | Pe| 7^ \Ph\- Finally, in Eq. ((10)). the integration over 
that part of T that lies in the potential-free region gives 
2el/v-F, where I represents the path length of the elec- 
tron travelling in the potential-free region. Putting all 
these pieces together, the total integral in Eq. (fTU|) can 
be expressed as S e — Sh = 2el/v-p, where 



I = I + l\\ sintf + 7T — — . (12) 



This result implies that for soft- wall billiards with a con- 
fining potential given by Eq. QJ, the Bohr-Sommerfeld 
approximation to the state counting function (Eq. ((5J) is 
still applicable if the classical geometric pathlength dis- 
tribution P(l) is modified to account for the potential 
contribution to the action (Eq. ((11)) ). The BS quantiza- 
tion therefore involves a modified pathlength distribution 
P(l). Accordingly, we perform a Monte Carlo simulation 
with typically 10 6 trajectories, where each trajectory is 
calculated by a forth-order Runge-Kutta integration. We 
determine its length I in the potential free region as well 
as the contribution from the displacement parallel to the 
soft wall Subsequently, I is transformed to I (Eq. (Q ) 
for each encounter with one of the soft walls. Note that 
the return probability P(l) obtained in this way differs 
from that calculated from the geometric length of the 
curved trajectories. The effective path length I (Fig. J3J)) 
associated with the action of the particle in the harmonic 
potential is longer than the geometric length l g . The rea- 
son for this surprising behavior is that while the individ- 
ual actions Sj_ and S± of the particle and the hole prop- 
agating in the parabolic potential are smaller than l g , 
with the geometric length l g and the constant momentum 
Pa. \v=o = the difference between the particle and the 
hole actions Sj_ — S± is larger than l g (p°l e — p± ). 
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FIG. 3: Comparison between the geometric path length dis- 
tribution P(l g ) (taking into account the length of the curved 
trajectories) and the modified path length distribution P(l) 
(Eq. (H2J). Both l g and I are given in units of L. The inset 
shows the geometry: Andreev billiard with one parabolic wall, 
W = 0.7L, a = 6.7, x = OAL. The letters (b)-(d) denote 
the lengths / of those bundles of classical orbits, at which the 
wave functions show enhancement in Fig. 0] 



III. NUMERICAL RESULTS 

In this section we compare the state counting function 
obtained from the exact numerical quantum calculations 
with that predicted by the BS approximation presented 
above. We consider a ballistic dot confined by (a) one or 
(b) three parabolically shaped potential walls as shown 
in Fig. Moreover, we present both particle as well 
as hole components of the wavefunction of eigenstates of 
these N-S hybrid systems. 



A. One soft wall 

We first consider an N-S system in which the normal 
dot is confined by one parabolic potential as shown in 
Fig. HJa). Numerically we have found that most classi- 
cal trajectories starting from the N-S interface will hit 
the parabolic wall once, before returning to the super- 
conductor (see e.g. the trajectory indicated in Fig. |5J). 
This feature allows us to test our semiclassical approach 
on a comparatively simple level before proceeding to the 
more complicated case of three parabolic walls. We find 
that the agreement between the exact quantum mechani- 
cal calculations for Nqm(s) and its semiclassical counter- 
part Nbs(z) is generally very good (see Fig. Ufa)), which 
proves that our semiclassical approach extended to An- 
dreev billiards with soft walls is well suited to describe 
these systems. A small discrepancy arises only at excita- 
tion energies above the cusp marked (b) in Fig. Ufa). We 
also note, that both curves show a distinct cusp struc- 
ture, which has also been found in other Andreev bil- 
liardsw^i^ These cusps can be understood from the 
BS approximation as being the consequence of the behav- 
ior of the path length distribution P(l) at certain path 



N(e) 



(a) 




FIG. 4: (color online) (a) Comparison between the quantum 
mechanical state counting function N(e) (solid red staircase) 
and the BS prediction, Nbs(s), (Eq. (0, dashed green line) 
for a N-S system shown in the upper left inset. The modu- 
lus square \un\ 2 and \vn\ 2 (see Eq. (0) of the three eigen- 
states marked in (a) are shown in (b-d) to illustrate the corre- 
spondence between continuous families of classical trajectories 
and the quantum mechanical wavefunctions. The parameters 
are M = 15, W = 0.7L, h = 0, a = 6.7, x 0.4L and 
A/E F = 0.02. 
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FIG. 5: (color online) (a) Quantum mechanical state counting 
function (red solid staircase) and BS approximation (green 
dashed line) for a = 5.5. The parameters are M = 30, W = 
L, a = 125, x = L and A/E F = 0.02. (b) - (d) show selected 
eigenstates (see text). 



lengths. For example, as shown in Fig. [31 P(l) has a peak 
at I = 2L which results in the most pronounced cusp of 
the state counting function at e/A = 0.63, as obtained 
from Eq. (0). This cusp, marked (b) in Fig. Ufa), cor- 
responds to the shortest classical orbits possible in the 
system. The electron leaves perpendicular to the N-S in- 
terface, is reflected back at the soft wall, and then reaches 
the N-S interface again. Such trajectories are similar to 
the stationary chords found in Ref. [l]I The correspond- 
ing peak in P(l) (also marked as (b) in Fig.|3J) represents 
a continuous family of trajectories (or bundles^) that 
feature the same topology and bouncing pattern and (al- 
most) the same length, but different initial conditions. In 
this case, the bundle is generated by a continuous change 
in the y coordinate at the N-S interface. The correspond- 
ing probability density of the wave function in Fig. 
displays a pronounced enhancement along the bundle of 
classical periodic orbits. A similar phenomenon has re- 
cently been found in pseudointegrable normal billiards^ 
and was called superscarring. 



According to Eq. ((SJ, for constant quantum number 
n, smaller excitation energies e/A correspond to longer 
classical orbits. This suggests that the excitation energies 
below the pronounced cusp at e/A = 0.63 correspond to 
classical orbits with lengths I > 2L. Indeed, for smaller 
e/A certain wavefunctions mirror bundles of longer clas- 
sical orbits. Some examples are shown in Fig.Qfc, d), (see 
also Fig. [31 for the correspondig path lenghts). The "su- 
perscarring" of the corresponding wavefunction near the 
bundles of short orbits also suggests that the BS approx- 
imation reflects the essential features of the underlying 
quantum dynamics. 



B. Three soft walls 

We now consider the more complicated case of N-S 
systems in which the normal dot is confined by three soft 
walls as shown in Fig. 

The agreement between the semiclassical BS approxi- 
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mation (Fig. E(a)) and the quantum mechanical calcula- 
tion for the state counting function is again very good. 
The position of the cusp in Nqm(s) can be predicted 
by the BS approximation to originate from the bundle 
of shortest trajectories. The corresponding peak in P(l) 
(shown in the lower right inset of Fig.Efa)) is located at 
I = 2.7 L resulting in a cusp at e/A = 0.49. The wave 
function of the energy level at the cusp (see Fig. Efb)) 
shows a density enhancement over a bundle of classical 
periodic orbits analogous to the one discussed for the dot 
with only one soft wall. As the shortest orbits do not ex- 
plore the upper and lower horizontal walls, this similarity 
is to be expected. 

Fig-EJcf) presents projections of an Andreev state that 
explores the additional soft walls. It mirrors in both the 
particle and hole components a bundle of classical trajec- 
tories all bouncing in sequence at all three soft walls. In 
line with the retracing approximation, their length I = 3 
is reflected in P(l) (marked as (d)). The "width" of the 
bundle (i.e. the continuous range of initial conditions) is, 
however, narrower than the bundle discussed in the pre- 
vious paragraph, resulting in a strongly reduced slope in 
Nbs(s). 

A very different scenario emerges for the Andreev 
eigenstate depicted in Fig. Efc). The particle and hole 
parts display a drastically different pattern pointing 
to the limitation of the ideal retracing approximation. 
While the hole part (x) is characterized by the "bounc- 
ing ball" type scar of trajectories bouncing at the up- 
per and lower horizontal wall, the electron part i^n(x) 
shows no distinct enhancements. Furthermore, integrat- 
ing |^n(x)| 2 and |^n(x)| 2 over the area of the normal 
dot we find that the probability of finding the quasipar- 
ticle in the hole state is 79% while the probability of the 
electron state is only 12% (and altogether 9% of finding 
the quasiparticle in the superconductor). Taking into 
account also that the hole component ^n(x) has a very 
low amplitude at the N-S interface, one can come to the 
conclusion that the hole part is only weakly coupled to 
the superconductor. The Andreev state of Fig. [5f c) can 
therefore be quantized by assuming that the particle is 
quasi-bound in the hole space where it spends most of its 
time with only infrequent and short excursions into the 
electron space, where it travels on trajectories scattering 
off the soft walls before returning to the hole space. This 
observation suggests that for certain individual Andreev 
states, semiclassical EBK quantization becomes possible, 
complementing the BS approximation for the smoothed 
state counting function. The almost vanishing density 
of |^n| 2 at the N-S interface allows us to assume that 
the wave function satisfies hard wall boundary conditions 
along the left wall of the cavity. Accordingly, the eigenen- 
ergies are to leading order determined by the quantized 
actions in the hole space in analogy to long-lived shape 
resonances. The EBK quantization conditions for the 



"isolated" hole state therefore read 



2p x W x 



2 Py Wy 



2^f2m*aE F 
V2ra*a£ F 



27rft 



( , V\ 



2^ (m+Hf) 



(13a) 
(13b) 



where W x and W y are the linear dimensions of the rect- 
angular region where V = 0, while p x and p y are the 
x and y components of the momentum of the particle. 
Here n and m are the quantum numbers characterizing 
the actions on the tori and \i\ and ji2 are the Maslov 
indices. The left hand sides of these two equations de- 
scribe the classical actions corresponding to the motion 
of the particle along the x (Eq. (|13a)0 and y (Eq. (|13b[0 
directions. 

The Maslov index of /ii = 3 in (Eq. (|13a|) ) originates 
from one reflection at the soft wall and one reflection 
at the hard wall which replaces the N-S interface, while 
ji2 = 2 in Eq. (|13b|) due to the reflections at the two soft 
walls. The two equations can be solved for p x and 
p y as functions of the two quantum numbers n, m. The 
eigenenergy of the N-S system in EBK quantization is 
then given by e = (p 2 -f p 2 )/(2ra*) — E^. By inspection 
of Fig. |5{c) we find for the hole part of the wave func- 
tion |^at| 2 three maxima along the x, and 35 along the y 
direction, respectively. Thus, the two quantum numbers 
in Eq. (T3J) are to be taken as n = 2 and m = 34, result- 
ing in | e/A | = 0.47. The excitation energy of the state 
shown in Fig. UJc) was found to be e/A = 0.448. This is 
a remarkable agreement, given that the mean level spac- 
ing 5n of the normal dot in units of A is 0.05, and taking 
into account that semiclassical methods usually cannot, 
to the leading order, resolve the energy spectra on much 
finer scale than the mean level spacing. The small differ- 
ence between the energies obtained from the BdG equa- 
tion and the semiclassical calculation can partly also be 
attributed to the fact that the wave function at the N- 
S interface is not exactly zero, i.e., a weak coupling of 
the normal cavity to the superconductor is present. This 
result implies that the Andreev state in Fig. [5f c) is, to 
a very good approximation, equivalent to an eigenstate 
of the hole in an isolated cavity. A similar effect of the 
decoupling of the wave function from the superconduc- 
tor has been observed in a system of a superconducting 
disk surrounded concentrically by a normal conductor™ 
In this hybrid system such states were called whispering 
gallery states. 

Finally, we investigate how a smaller a, correspond- 
ing to a shallower soft wall, affects the agreement with 
our semiclassical description. We consider the same sys- 
tem as in Fig. [51 only the parabolic slope a (or steepness 
of the wall) is reduced from 125 to 5.5. In such a sys- 
tem, the particle spends more time in the region where 
V 7^ 0. The resulting state counting function is shown 
in Fig. EJa). The cusp structure is predicted very accu- 
rately by our BS approximation. The most pronounced 
cusp, marked (c) in Fig. |B{a), is located at an energy 
of e/A = 0.37. The position of this cusp can again be 
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FIG. 6: (color online) (a) Quantum mechanical state counting 
function (red solid staircase) and BS approximation (green 
dashed line) for the same system as in Fig. |5] except for a 
smaller value of a = 5.5 (curvature of parabolic wall), (b) — (d) 
show selected eigenstates (see text). 

determined by the peak in the path length distribution 
P(l) at I = 3.3L, corresponding to the bundle of short- 
est trajectories. Note, however, that the non-geometric 
shift in length introduced by Eq. (|12[) is 1.3L, and is now 
comparable to the system size. 
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